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Abstract 

Given a chaotic dynamical system and a time interval in which some quantity takes an 
unusually large average value, what can we say of the trajectory that yields this deviation? 
As an example, we study the trajectories of the archetypical chaotic system, the baker's map. 
We show that, out of all irregular trajectories, a large-deviation requirement selects (isolated) 
orbits that are periodic or quasiperiodic. We discuss what the relevance of this calculation may 
be for dynamical systems and for glasses. 



1 Introduction 



In this paper we discuss, with the aid of a simple example, how order arises from chaos when we 
select the solutions that minimize some quantity. The motivation for this work is twofold, one 
related to dynamical systems, and the other to glasses. 

Consider first a chaotic dynamical system, for example a liquid flowing past an obstacle, at large 
Reynolds numbers. The wake behind the obstacle is chaotic, and the pressure it creates may exhibit 
considerable fluctuations. We may ask the question: given that during a certain time-interval the 
average pressure has been abnormally large, what can we say about the velocity field during that 
time? Is it reasonable to assume that, if the high pressure is to last long, then the fluid phase- 
space trajectory has to have during that time some regularity that is absent in a typical chaotic 
trajectory? What we are asking, then, is a question on the large deviation functions of a chaotic 
system: are the trajectories that sustain extreme situations special, and in a sense simpler? 

A second motivation is related to the question of the nature of the order in an ideal, equilibrium, 
glass state - if, of course, such a thing exists. An approach to the glass problem [1] is to consider 
a free-energy density functional in (i-dimensional space [2] : 

F[p(x)] = J d d x p[lnp(x) ~l]- l -Jd d * dV [p(x) - Po ]C(x - x')[p(x') - Po ] (1) 

where C(x — x', p a ) is the liquid direct correlation function at average density p Q . We look for the 
'local' free energy minima that satisfy: 



6F[p(x)] 



= lnp(x) - J C(x - x', Po )[p(x') - Po ] = (2) 



At low average densities p a , the spatially constant 'liquid' solution dominates. As the density 
increases, a periodic, 'crystalline' solution appears. What is interesting from the glassy point 
of view, is that in the high density regime, there appear also many non-periodic 'amorphous' 
solutions. Each one of these is supposed to represent a metastablc glassy state. These states are 
local minima of Q satisfying Q. Beyond the mean-field approximation, high free-energy solutions 
are unstabilized by nucleation; if we wish to model the realistic situation we should look for the 
solutions that are deepest in free energy - excluding, of course, the crystalline state. We have then 
a situation analogous to the one described above for chaotic systems: of all the solutions satisfying 
the " equations of motion" ([2| , we have to look for the configurations that minimize globally the free 
energy (JT|), and ask whether these "large deviation" solutions have special regularity properties. 

This analogy between dynamical systems that are chaotic in time, and glassy systems that are 
chaotic in space, was pointed out many years ago by Ruelle |3J. It has also appeared in the theory 
of charge-density waves |U[S], in particular in the Frenkel-Kontorova model. It turns out that the 
local energy minima of the model are given by the trajectories of the "standard map", which has 
both regular and chaotic orbits. However, when one restricts to the lowest energy minima, the 
trajectories selected are quasiperiodic. 

In a recent paper it was argued [5] that ideal glasses should exhibit a form of geometrical order 
that is more general than periodic and quasiperiodic, characterized by the fact that local motifs 
of all sizes are repeated often. For a system that is not completely ordered, this assumption leads 
to the definition of a correlation length, defined as the largest size of patterns that repeat more 
frequently than they would in a random case. The claim is then that this correlation length should 
diverge in an ideal glass. But frequent repetitions of patterns would not occur for a generic solution 
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of a spatially chaotic glass. If order in the solution is to happen at all it is through a mechanism 
as described above: choosing configurations that globally minimize energy would then force the 
system to favor some patterns that would consequently be repeated often. 

In the present work we treat a simple example. The system studied has the advantage that it is 
purely chaotic, we know that there are no islands of regular dynamics in phase-space. If the large- 
deviation bias selects orbits which have some regularity, it does so in the most unfavorable situation 
in which these orbits are isolated and unstable. The organization of this paper is as follows: In 
section [2j we present the baker's map and the functional used in our model. In section [3j the 
trajectories that extremize the functional are found and their order properties are discussed. 



2 The model 

In this section we briefly review the baker's map and some of its properties. We then define the 
functional used in our model and whose minimal trajectories will be studied later. 

2.1 The baker's map 

The baker's map is a simple two-dimensional chaotic map acting on the unit square. The dynamical 
system with discretized time t, associated with the baker's map, is defined for q t G [0, 1] and 
Pt € [0,1] by the area-preserving equations 



Qt+i 
Pt+i 



2q t 
1 
2 



{Pt+ [2q t \), 



(3) 
(4) 



where [2qt\ is the integer part of 2q t . As shown by figure [T] the phase square is stretched in the q 
direction and squeezed it in the p direction; the baker then cuts the right part and puts it on top 
of the left one. 














Figure 1: One iteration of the baker's map. 

The baker's map can also be understood as a shift operator. Let us consider the binary rep- 
resentations qt = O.abc . . . and p t = O.uvw . . . where a, b, c, u, v, w, . . . are either or 1. If the 
position (qt,Pt) of the system at time t is written as the string "... wvu-abc. . .", the state of system 
at time t + 1 is obtained by shifting the central dot to the right by one position, that is, the string 
becomes ". . . wvua ■ be . . . " and we have qt = O.bcd . . . and p t — O.auv .... 

The baker's map is deterministic: if the system is started at initial conditions (<?0;Po)j equa- 
tions pi) and (H deter mine its position at all future (and past) times. However, the system is 
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chaotic in the sense that a small uncertainty in the initial conditions grows exponentially in time 
(the shift operator doubles it at each iteration); when the uncertainty reaches the size of the phase 
square (here, this size is unity), the position of the system becomes unpredictable. 



l 








Figure 2: Representation of the periodic orbit {ri4, . . . , no} = {0, 0, 1, 0, 1} in the phase square 
(left). The arrow shows the direction of the trajectory. An irregular orbit that never is 
an extremal trajectory, for any g* (right). 

Although almost all trajectories are aperiodic, there are many initial conditions leading to peri- 
odic orbits. Considering the shift operator interpretation and remembering that the binary repre- 
sentation of a rational number has a periodic pattern, one sees that the periodic orbits are related to 
rational numbers. In figure [2] (left), the rational number 5/31 ( whose binary representation is the 
periodic string 0. 00101 00101 00101 . . . ) was chosen as the initial position go (the top-left point). 
The ordinate po is such that its binary representation is the mirror of go's binary representation. 
After five units of time, the dot shifts five times to the right and one has gs = go; thus, the orbit 
is periodic with period 5. A periodic orbit can be described by the binary digits of the periodic 
pattern in the binary representation. In our example, the orbit is {714,..., no} = {0,0,1,0,1}. 
Obviously, a cyclic permutation of these five digits doesn't change the orbit. 

Strictly speaking, the binary representation of a rational is not always periodic, but the non- 
periodic part is located near the dot and has a finite length. For example, the number go = 5/31+1/4 
is written as 0. 011 01001 01001 . . . After many iterations of the map, the system forgets the 1/4 
term in g (its second most significant bit), because it is shifted to less and less significant places in 
the binary representation of p t ; the initial po is forgotten as well. It is clear that as long as the initial 
position go is rational, the trajectory falls into a periodic orbit after some time. If a go is chosen 
randomly and uniformly in the unit interval, it is irrational. Indeed, although rational numbers are 
dense in the unit interval, they form a set with zero measure. As a consequence, typical trajectories 
are not periodic. Because the baker's map is ergodic, such a typical trajectory fills up uniformly the 
phase square. Furthermore, a small perturbation of an initial condition associated with a periodic 
orbit yields a chaotic orbit. 

The baker's map has two unstable fixed points which correspond to periodic orbits with period 
one. These two orbits are {0} and {1}, respectively the points located at the origin (0,0) and the 
opposite corner (1,1). 
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2.2 Large deviations 



Let us introduce in our model an observable which is a functional of the trajectories. Our choice is 
arbitrary, we are motivated here by simplicity. We consider a functional parametrized by a number 
g, in [0,1]: 

1 T_1 

h[{qt},q*]= lim ^ $> t - g*) 2 - (5) 

t=0 

The observable h can be interpreted geometrically as the average squared distance over time between 
the trajectory and a reference vertical line located at g*, as schematically represented in figure [2] 
For each value of the parameter q„, there is one particular trajectory that minimizes the average 
distance h. Note that h is invariant by the transformation (q t , g*) — ► (1 — qt, 1 — g*) so if {q t } is an 
extremal trajectory for g* < i, {1 — qt} is an extremal trajectory for 1 — g* > 1. Hence, we can 
restrict the study to < g* < i. Although the baker's map is bidimensional, the problem is actually 
one-dimensional only, because the dynamics on the q axis are not coupled to the one on the p axis, 
as shown by the equation ^ . The two-dimensional generalization of the observable ^ would be 
the average squared distance between the trajectory and a reference point located at (g*,p*): 

lim ±J2[(q t -q*) 2 + (p t ~p*) 2 ]. (6) 

1 — >oo ± — 

t=0 

The extremal trajectories for (g*,p*) of the observable ^ are easily deduced from the extremal 
trajectories of the one-dimensional observable ^ for the parameter \(q* + p*). It will indeed be 
clear later that the extremal trajectories are left invariant by the transformation that swaps the 
q and p axis; this transformation is equivalent to reversing the time. The orbit in figure [2] is for 
example the minimal trajectory for the chosen value of q* (and actually for any g* in the whole 
range [0.385,0.411]) and has the (p <-> q) symmetry, as do the orbits of figure [6j 

The worse state in the sense of minimization of h is one of the two fixed points. If g* < 1, this 
maximal trajectory is the fixed point (1,1) and we have h — (1 — g*) 2 . On the other hand, the 
distribution of the qt for a chaotic trajectory is uniform in the unit interval and the mean distance 
is calculated by the continuous sum 

f 1 1 

frchaotic(<?*) = J (g- g*) 2 dg = - - g*(l - g„). (7) 

Figure [3] shows the energies of the minimal orbit (which we shall determine below) , the maximal 
orbit and a typical (i.e. chaotic) orbit. It is worth noting that chaotic orbits never minimize or 
maximize the functional. In the next sections, the minimal trajectories of h are made explicit and 
found to be periodic for most values of g* and quasi-periodic for some (zero measure) set of values 
of g*. 



3 Large-deviation function as a free energy of a ID lattice 
gas model 

3.1 Mapping to the particle model 

In this section, we will show that the functional ^ can be interpreted as the Hamiltonian of particles 
on a one-dimensional lattice with single occupancy and a translationally invariant interaction. That 
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Figure 3: Value of the observable ^ for the minimal trajectory (lowest curve), the maximal 
trajectory (highest curve) and a typical trajectory (middle curve). 



is, 



h = — ^2 U\i-j\ n i n j - VP + h o 



where p is the particle density, and \i acts as a chemical potential. The interaction will be shown 
to decay exponentially with distance. 

As stated earlier, a point on a periodic orbit (of period N) is represented by the binary sequence 
0.njv-i-"^i^o^Ar— I--- = {npf— i, . . . , no}. Note that this number may be written as 

{njv-i, • • • , n } = N l _ (2 N ~ 1 n N - 1 H h 2 1 n a + n Q ) . 



G 



We define a as the cyclic left shift operator (or left cyclic permutation) to be applied on the 
digits rii of the period, that is 

cr{riN-i, ■ ■ ■ ,ni,n } = {n N ^ 2 , ■ ■ ■ ,n ,n N -i}. 

The i th iteration of a may be written as 

N-l 



3=0 

Here and in what follows [•] denotes the exponents and indices modulo N. 

Taking into account the periodicity and expanding, the observable ^ now reads 

h = u-fip+h , (8) 

where 

1 

P= T7 E M = 2?*, h = qt- (9) 



and with u given by 



i jk 



njny +r ] (10) 



In the above, all the indices are summed from to N — 1. Shifting the exponents on the right hand 
side as i — > [i — j] we get: 

« = ^ E ^ ^ = T^rp E 2[4]2[r+l] en) 

jr i 

In the Appendix we show that 
3.2 The long-orbit limit 

Consider now a very long chain of size M, particles with occupation numbers m, and an interaction 
given by 

1 M i M 1 M 

^=^E U\i-j\ n i n i ~ Yl^ nt + ha ^ Ml E 2 H *~ j| -2g*p+ (q*) 2 (13) 

i,j=l i=l ' i,j=l 

that is, [/ r = |2~' r '. Let us compute the energy of a periodic configuration of period N << M. 
We can re-express the unrestricted sum (13) as a sum over the first period as in Figure |4j adding 
first each position r within the first period with its repetitions in all other periods. We may express 
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these contributions as two geometric series, corresponding to the particles to the right and to the 
left, respectively: 



1 

3 

1 1 



3(1 



| 2~ r + 2~ (JV+r) + 2~ ( - 2N+r) + ..." + ~2~ {N - r) + 2-< 2JV - r ) + 2-( 3Ar - r ) + ..." | 

i oTV — r i or 

[2- r + 2- f 1= Vr = ^ (14) 



)-N 



Thus, (11 1 is recovered. This means that we may consider, without loss of generality, a very long 
periodic chain of size M, and the energy of periodic configurations of periods N << M will be 
correctly computed by the interaction (13 1. This, in turn means that our large functional (p| can 



be computed for large times on the basis of the particle model (13 1 





N 










► 














r , N+r, N+2r, ... 


N-r, 2N-r, 3N-r, 





Figure 4: Counting the energy of a periodic configuration as the sum of two geometric series, 
corresponding to the set of periods to the right, and to the left, respectively. 



3.3 Large deviation function 

Given a time T, we express the probability of observing a value of the quantity defined in (|5| as: 

P(h{{q uq *} = h) = e- T ^ ; (15) 

this defines £(/i). Alternatively, we may express the same information in the large deviation function 
/(A) : 

e- Tf(X) = J d\ P(h[{q t , 9 »] = h) e~ AT h (16) 

For long times T, we have that £ and / are related by a Legendre transform. In fact, we recognize 
/(A) as the free energy density of the particle model at inverse temperature A. In particular, the 
ground state (corresponding to infinite A) is given by the trajectory minimizing h. 



8 



3.4 The ground state of the lattice particle model 

Let us start by considering the minimum of the energy ( |13[ ) : 



min{/i} = min 



5> 



r] -2q,p+{q,f 



.)>' 



mm 

p 



{U min (p)-2q*p+(q*) 2 } (17) 



where U mm (p) is the smallest energy for a fixed configuration density. Although it is tempting to 
determine the value of p of the extremal configuration by 



dU min (p) 
dp 



2q* 



(18) 



this must be treated with care, because as we shall see, the derivative is discontinuous. Note again 
that q* plays the role of a chemical potential. 

Our job is now to determine U mm (p). This problem was solved years ago [TJ [H] for long-range 
(U r > Vr) interactions with U r convex downwards, as is the case here. 

For the moment, let us forget the lattice. If the particles can take any position on a one- 
dimensional line with periodic boundary conditions, the ground state is simply a crystal, that is, 
each particle is separated from its neighbour by the same distance d = 1 /p. The condition of 
convexity of the interactions can be understood as a stability condition. Indeed, if one considers a 
first-neighbor concave interaction as shown in Figure [5] (left), a displacement of one particle in the 
crystal would decrease the energy of the system: the system is not stable. On the contrary, if the 
first-neighbor interaction is convex as shown on the figure [5] (right), a displacement of one particle 
would increase the energy of the system. 

A 

• e-« • 



d + 8 d 



Au<0 




d + S 




Au>0 



d + S 



Figure 5: Change Am of the energy of a cristal configuration after the displacement of one particle 
for a concave interaction (left) and a convex interaction (right). 

For particles constrained to lie on a lattice, if p is such that the mean distance 1/p is an integer, 
then the ground state is again periodic with particles separated by this same distance l/p. In the 
baker's map representation, the trajectory corresponding to such a crystal is very regular. Figure [6] 



9 



shows the p — |, p = | and p = | minimal trajectories. As 1/p becomes larger, the trajectory 
spends more and more time near the fixed point (0, 0). 

If l/p is not an integer, then it lies between the two successive integers d = [l/p\ and d+1, and 
the distance between any two consecutive particles will be either d or d+1. There is a competition 
between the solution with density 1/d and the one with density l/(d+l). This is seen in the baker's 
map, where Figure [6] shows that the trajectory for p = 5/17 spends some time in the p = 1/4 (period 
4) solution and some time in the p = 1/3 (period 3) solution. The competition between the two 
solutions 1/d and l/(d+l) is also seen in Figure [7j In this figure, the 1/d solutions are seen to 
dominate for some large range of values of q* . When varying the parameter continuously in order 
to go from the 1/d solution to the l/(d+ 1) solution, the period takes many intermediate values 
"interpolating" between periods 3 and 4. 

Let us now turn to the full solution, which was originally given in terms of continued fractions 
[8] and admits a graphical construction related to that employed in quasicrystals[9j [10] which we 
will describe. Consider the plane and the lattice of points with integer coordinates, as shown in 
figure [8] and the line with slope p starting at the origin (0, 0). The algorithm to solve the problem is 
quite simple: for each integer abscissa k, compute the ideal ordinate kp and round it to the closest 
smaller integer y^, then, mark the point (k,yk). The ground state corresponding to density p is 
simply given by the sequence of binary digits — yk — J/fc-i, that is, there is a particle in site k 
in places where the black pixels have a step upwards. In fact, the line can start at any arbitrary 
point; this leads to different sequences with the same energy. These sequences have the form 

n k =[kp + (j ) \-[(k-l)p + (f>\ , (19) 

where is a phase factor related to the y-intercept of the line. Figure [8] shows the determination 
of the ground state for p = 3/7. One can read directly the ground state by looking at the ordinates 
of the black pixels. 

Following Aubry [TT], one can obtain an exact expression for the energy U mm (p). We need 
to introduce first another useful representation of the ground state. If Xj is the position the j-th 
particle, we have 

Xj =[- + const\ . (20) 

Next, we rewrite the interaction term as 

U = J2 u r (21) 

r 

where the u r is the average contribution of all pairs of r th neighboring particles 

u r = p(U(xj +r — Xj)). (22) 

and (•) denotes averaging over all such pairs, in the limit M — > oo. Let us derive an explicit 
expression for the u r . The distance between two particles separated by r — 1 others, is 

cither L r /° _1 J or + 1- If n r is the fraction of pairs having the value [ r P _1 J > an d (1 — 7Tr) of 

having [rp~ J+l, (xj +r — xj) is by definition ir r \rp~ 1 \ + (1 — 7r r ) (|_?"P _1 J + !)■ Using the fact that 

{x j+r - Xj) = (x j+r - (j + r)/?" 1 ) - (xj - jp- 1 ) + rp- 1 = rp^ 1 , 

we conclude that 

7r r = 1 - (rp- 1 - [rp- 1 ^) . 
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Figure 6: Examples of minimal orbits for different values of q« . The density of l's p of the orbits 
are indicated. 
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Figure 7: The period of the the minimal trajectory as a function of g*. The fractions specify 
the density p of particles in the corresponding configuration. Only the largest steps are 
indicated. 
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Figure 8: Determination of the ground state of the lattice Hamiltonian for p = 3/7. The arrows 
indicate the position of the particles. 
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We thus obtain an explicit expression of the interaction term: 



U min = pJ2 [irrUdrp- 1 ]) + (1 - n r ) U{\rp- l \ +1)]^E( 2 ' - L^J)) pU{\rp- l \). 

r r 

(23) 

jjmin - g con tinuous everywhere and its derivative 

u'(p) = = J2 (1 + IL^-'J) Udrp- 1 }) (24) 

r r 

is an increasing function, as shown in Figure [9] 

The density can be determined from the knowledge of following equation (18 1: the y-axis on 
figure [9] is 2g*. When the level 2q* hits on an interval between two values, the lower step is chosen. 
The density p is thus given by Figure |10| Once p is known, the configuration can be obtained either 
graphically from Figure [8j or with equation (20 1. The value of the 'energy' can be obtained from 
Equations (17 1 and p3| ), or the first of figures|9 

If 2g* falls between two steps of the u'(p) curve, i.e. there exists a p such that u'(p~) < p < 
u'(p + ), then p is rational and is the density of the ground state; the configuration of particles is 
consequently periodic. The width of the interval of corresponding to this density is given by the 
discontinuity of the energy's slope 25q Sf (p) — u'(p + ) — u'(p~) If, on the contrary, there exists a p 
such that u'(p) = 2q*, then p is irrational: the ground state has density p and the configuration of 
particles is quasi-periodic. 

There exists a smallest integer r such that rp^ 1 = s where s is an integer. For this r and 
for all its multiples kr, u' kr (p) is discontinuous and J2 k r (u' kr (p + ) — u' kr (p~)) is the sum of the 
discontinuities. We obtain 



00 

H p = 1) = 3 B fcs )2- fes 



sJ 3^ v ' 3 (2 s -l) 2 ' 

k=l y ' 

Note that the result does not depend on the numerator r. 

Although a quasi-periodic ground state can be obtained by choosing the right , the set of such 
particular values has zero measure. Indeed, the sum of all discontinuities of u'(p) for < p < 1 
rational, is 

00 00 

a = y, - u '(<°")) = 3 E E ^)(M2- fcs , 

p=§e[o,i] «=i k=i 

where <^(s) is the Euler function equal to the number of integers in [1, s — 1] that are coprime to s. 
Using the change of variable u — ks and a known property of the Euler function, we obtain 

1 00 -.00 . 

^e^e^He^H- 

«=1 k\u u=l 

We have proved that A = | = u'(l) — it'(0), that is, the probability that a fixed q* falls between 
steps is one. In other words, the probability that the density is rational is one and the set of values 
of associated with quasi-periodic orbits has zero measure. 



The solution p(q*) plotted in figure 10 has a devil's staircase structure: it is a continuous curve, 



constant on some interval when p is rational. 



13 



0.6 



0.5 



0.4 



■a. 0.3 



0.2 



0.1 




Figure 10: The density p of the ground state plotted versus the chemical potential /i = 2q*. This 
curve is a complete devil's staircase. 
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4 Concluding remarks 



We have studied the trajectories of the baker's map that extremize the average squared distance 
to a line q = q*. The minimal trajectory is periodic for most values of g*; and quasi-periodic for 
a zero-measure set of values of q*. For all values of q*, the minimal trajectory exhibits more order 
than the generic chaotic trajectories. This does not mean that the trajectories involved are not 
chaotic in the sense of being insensitive to boundary conditions: a slight perturbation of the initial 
condition of an extremal trajectory leads to a new trajectory that departs away from the periodic 
one, and ultimately explores ergodically all phase-space - thus spoiling its extremal properties. 

The ordered structure of the minimizing state highly depends on fact that the corresponding 
Hamiltonian in the ID lattice gas has infinite range. In the case in which the Hamiltonian of the 
ID lattice gas has a finite range ( U r = for r > R - one may think about a hard-sphere model) , 
for small enough density any configuration with interparticle distance > R minimizes the energy; 
and the vast majority of configurations is not ordered. This is easy to understand for the baker's 
map: such an interaction corresponds to minimizing a functional that does not depend on all digits 
of the positions, i.e. it is a functional in which the trajectories contribute in a degenerate manner. 

In the other extreme, functionals that are long-range in time may also yield disordered extremal 
trajectories. For the baker's map, the observable to be minimized could be 



where [2q t \ is the t-th binary digit of go- This is a long-range antifcrromagnet: the sum thus has 
disordered, degenerate ground states. Note that this observable is not local in the sense that the 
position at time t is coupled to the position at different time t'. 

In conclusion, we have presented a simple example illustrating that those trajectories of a map 
which extremize some functional may exhibit structure absent in the entire ensemble of trajectories. 
This is reminiscent of the Frankel-Kontorova model, where the ground states possess quasiperiodic 
translational orderjH while other stationary solutions do not. It is tempting to conjecture that 
such behavior is the rule rather than the exception, and may therefore be relevant to systems which 
are chaotic in cither the spatial or temporal domain. 
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Appendix 

Let us start with 



We first determine Uq N ^: 



^^T^T^E 2 ^ ( 25 ) 



12 N + l 

'o — 



u ° = nw=i- (26) 



Next, let us express U^ + [ in terms of Ur . We split the sums in 



N-l 



U r+ i 



^^E 2 ' 2 ^ 11 ( 27 ) 



iV-l N-r-2 N-l 



i=0 

as 

E = E + E = E + E • ( 2 8) 

i=0 i+r+KJV i+r+l>N i=0 i=iV— f — 1 

so that 

/N-r-2 N-l \ 

^ = 72^^ E 2*2^ + £ 2*2^^ (29) 
(here the exponents are not necessarily modulo N). The first and second sums read: 

N-r-2 /N-r-1 \ 

g 2 l 2 J+r+1 =21 ^ 2 l 2 J+r+1 - 2 N - r - 1 2 N (30) 



i=Q \ i=Q / 

and 

N-l / iV-l \ 

jjV-r-1 



N-l / iV-l \ 

T2 l+r+1 - N =21 ^ 2 t 2 l+r - Ar 

i—N—i — 1 \i=N-r / 



or: 



with C = 1/(1 - 2 _JV ). Putting 
the recurrence equation becomes 



(31) 



so that equation ( 29 ) can be written 

C/W = 2C/r ) - (2] ^ lp (2^- 1 2^-2^- 1 ) (32) 

oJV-i — 1 

ulll = 2 ^ N) ~w^T (33) 



(2 r+1 U^l)=A(2 r U^)-C (34) 



6 r = 2 r uW - £ (35) 



6 r+ i = Ab r (36) 
17 



which implies that b r — &o4 r , with bo — Uq — leading to 
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